High performance Legionella pneumophila source attribution using genomics-based machine learning classification.

Fundamental to effective Legionnaires' disease outbreak control is the ability to rapidly identify the environmental source(s) of the causative agent, Legionella pneumophila. Genomics has revolutionized pathogen surveillance, but L. pneumophila has a complex ecology and population structure that can limit source inference based on standard core genome phylogenetics. Here, we present a powerful machine learning approach that assigns the geographical source of Legionnaires' disease outbreaks more accurately than current core genome comparisons. Models were developed upon 534 L. pneumophila genome sequences, including 149 genomes linked to 20 previously reported Legionnaires' disease outbreaks through detailed case investigations. Our classification models were developed in a cross-validation framework using only environmental L. pneumophila genomes. Assignments of clinical isolate geographic origins demonstrated high predictive sensitivity and specificity of the models, with no false positives or false negatives for 13 out of 20 outbreak groups, despite the presence of within-outbreak polyclonal population structure. Analysis of the same 534-genome panel with a conventional phylogenomic tree and a core genome multi-locus sequence type allelic distance-based classification approach revealed that our machine learning method had the highest overall classification performance-agreement with epidemiological information. Our multivariate statistical learning approach maximizes the use of genomic variation data and is thus well-suited for supporting Legionnaires' disease outbreak investigations.IMPORTANCEIdentifying the sources of Legionnaires' disease outbreaks is crucial for effective control. Current genomic methods, while useful, often fall short due to the complex ecology and population structure of Legionella pneumophila, the causative agent. Our study introduces a high-performing machine learning approach for more accurate geographical source attribution of Legionnaires' disease outbreaks. Developed using cross-validation on environmental L. pneumophila genomes, our models demonstrate excellent predictive sensitivity and specificity. Importantly, this new approach outperforms traditional methods like phylogenomic trees and core genome multi-locus sequence typing, proving more efficient at leveraging genomic variation data to infer outbreak sources. Our machine learning algorithms, harnessing both core and accessory genomic variation, offer significant promise in public health settings. By enabling rapid and precise source identification in Legionnaires' disease outbreaks, such approaches have the potential to expedite intervention efforts and curtail disease transmission.

High-performance enrichment-based genome sequencing to support the investigation of hepatitis A virus outbreaks.

IMPORTANCE: This proof-of-concept study introduces a hybrid capture oligo panel for whole-genome sequencing of all six human pathogenic hepatitis A virus (HAV) subgenotypes, exhibiting a higher sensitivity than some conventional genotyping assays. The ability of hybrid capture to enrich multiple targets allows for a single, streamlined workflow, thus facilitating the potential harmonization of molecular surveillance of HAV with other enteric viruses. Even challenging sample matrices can be accommodated, making them suitable for broad implementation in clinical and public health laboratories. This innovative approach has significant implications for enhancing multijurisdictional outbreak investigations as well as our understanding of the global diversity and transmission dynamics of HAV.

Genomic surveillance for antimicrobial resistance - a One Health perspective.

Antimicrobial resistance (AMR) - the ability of microorganisms to adapt and survive under diverse chemical selection pressures - is influenced by complex interactions between humans, companion and food-producing animals, wildlife, insects and the environment. To understand and manage the threat posed to health (human, animal, plant and environmental) and security (food and water security and biosecurity), a multifaceted 'One Health' approach to AMR surveillance is required. Genomic technologies have enabled monitoring of the mobilization, persistence and abundance of AMR genes and mutations within and between microbial populations. Their adoption has also allowed source-tracing of AMR pathogens and modelling of AMR evolution and transmission. Here, we highlight recent advances in genomic AMR surveillance and the relative strengths of different technologies for AMR surveillance and research. We showcase recent insights derived from One Health genomic surveillance and consider the challenges to broader adoption both in developed and in lower- and middle-income countries.

Bioinformatic investigation of discordant sequence data for SARS-CoV-2: insights for robust genomic analysis during pandemic surveillance.

The COVID-19 pandemic has necessitated the rapid development and implementation of whole-genome sequencing (WGS) and bioinformatic methods for managing the pandemic. However, variability in methods and capabilities between laboratories has posed challenges in ensuring data accuracy. A national working group comprising 18 laboratory scientists and bioinformaticians from Australia and New Zealand was formed to improve data concordance across public health laboratories (PHLs). One effort, presented in this study, sought to understand the impact of the methodology on consensus genome concordance and interpretation. SARS-CoV-2 WGS proficiency testing programme (PTP) data were retrospectively obtained from the 2021 Royal College of Pathologists of Australasia Quality Assurance Programmes (RCPAQAP), which included 11 participating Australian laboratories. The submitted consensus genomes and reads from eight contrived specimens were investigated, focusing on discordant sequence data and findings were presented to the working group to inform best practices. Despite using a variety of laboratory and bioinformatic methods for SARS-CoV-2 WGS, participants largely produced concordant genomes. Two participants returned five discordant sites in a high-Cτ replicate, which could be resolved with reasonable bioinformatic quality thresholds. We noted ten discrepancies in genome assessment that arose from nucleotide heterogeneity at three different sites in three cell-culture-derived control specimens. While these sites were ultimately accurate after considering the participants' bioinformatic parameters, it presented an interesting challenge for developing standards to account for intrahost single nucleotide variation (iSNV). Observed differences had little to no impact on key surveillance metrics, lineage assignment and phylogenetic clustering, while genome coverage <90 % affected both. We recommend PHLs bioinformatically generate two consensus genomes with and without ambiguity thresholds for quality control and downstream analysis, respectively, and adhere to a minimum 90 % genome coverage threshold for inclusion in surveillance interpretations. We also suggest additional PTP assessment criteria, including primer efficiency, detection of iSNVs and minimum genome coverage of 90 %. This study underscores the importance of multidisciplinary national working groups in informing guidelines in real time for bioinformatic quality acceptance criteria. It demonstrates the potential for enhancing public health responses through improved data concordance and quality control in SARS-CoV-2 genomic analysis during pandemic surveillance.

The two-component system WalKR provides an essential link between cell wall homeostasis and DNA replication in Staphylococcus aureus.

Among the 16 two-component systems in the opportunistic human pathogen Staphylococcus aureus, only WalKR is essential. Like the orthologous systems in other Bacillota, S. aureus WalKR controls autolysins involved in peptidoglycan remodeling and is therefore intimately involved in cell division. However, despite the importance of WalKR in S. aureus, the basis for its essentiality is not understood and the regulon is poorly defined. Here, we defined a consensus WalR DNA-binding motif and the direct WalKR regulon by using functional genomics, including chromatin immunoprecipitation sequencing, with a panel of isogenic walKR mutants that had a spectrum of altered activities. Consistent with prior findings, the direct regulon includes multiple autolysin genes. However, this work also revealed that WalR directly regulates at least five essential genes involved in lipoteichoic acid synthesis (ltaS): translation (rplK), DNA compaction (hup), initiation of DNA replication (dnaA, hup) and purine nucleotide metabolism (prs). Thus, WalKR in S. aureus serves as a polyfunctional regulator that contributes to fundamental control over critical cell processes by coordinately linking cell wall homeostasis with purine biosynthesis, protein biosynthesis, and DNA replication. Our findings further address the essentiality of this locus and highlight the importance of WalKR as a bona fide target for novel anti-staphylococcal therapeutics. IMPORTANCE The opportunistic human pathogen Staphylococcus aureus uses an array of protein sensing systems called two-component systems (TCS) to sense environmental signals and adapt its physiology in response by regulating different genes. This sensory network is key to S. aureus versatility and success as a pathogen. Here, we reveal for the first time the full extent of the regulatory network of WalKR, the only staphylococcal TCS that is indispensable for survival under laboratory conditions. We found that WalKR is a master regulator of cell growth, coordinating the expression of genes from multiple, fundamental S. aureus cellular processes, including those involved in maintaining cell wall metabolism, protein biosynthesis, nucleotide metabolism, and the initiation of DNA replication.

A statistical genomics framework to trace bacterial genomic predictors of clinical outcomes in Staphylococcus aureus bacteremia.

Outcomes of severe bacterial infections are determined by the interplay between host, pathogen, and treatments. While human genomics has provided insights into host factors impacting Staphylococcus aureus infections, comparatively little is known about S. aureus genotypes and disease severity. Building on the hypothesis that bacterial pathoadaptation is a key outcome driver, we developed a genome-wide association study (GWAS) framework to identify adaptive mutations associated with treatment failure and mortality in S. aureus bacteremia (1,358 episodes). Our research highlights the potential of vancomycin-selected mutations and vancomycin minimum inhibitory concentration (MIC) as key explanatory variables to predict infection severity. The contribution of bacterial variation was much lower for clinical outcomes (heritability <5%); however, GWASs allowed us to identify additional, MIC-independent candidate pathogenesis loci. Using supervised machine learning, we were able to quantify the predictive potential of these adaptive signatures. Our statistical genomics framework provides a powerful means to capture adaptive mutations impacting severe bacterial infections.

The importance of utilizing travel history metadata for informative phylogeographical inferences: a case study of early SARS-CoV-2 introductions into Australia.

Inferring the spatiotemporal spread of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) via Bayesian phylogeography has been complicated by the overwhelming sampling bias present in the global genomic dataset. Previous work has demonstrated the utility of metadata in addressing this bias. Specifically, the inclusion of recent travel history of SARS-CoV-2-positive individuals into extended phylogeographical models has demonstrated increased accuracy of estimates, along with proposing alternative hypotheses that were not apparent using only genomic and geographical data. However, as the availability of comprehensive epidemiological metadata is limited, many of the current estimates rely on sequence data and basic metadata (i.e. sample date and location). As the bias within the SARS-CoV-2 sequence dataset is extensive, the degree to which we can rely on results drawn from standard phylogeographical models (i.e. discrete trait analysis) that lack integrated metadata is of great concern. This is particularly important when estimates influence and inform public health policy. We compared results generated from the same dataset, using two discrete phylogeographical models: one including travel history metadata and one without. We utilized sequences from Victoria, Australia, in this case study for two unique properties. Firstly, the high proportion of cases sequenced throughout 2020 within Victoria and the rest of Australia. Secondly, individual travel history was collected from returning travellers in Victoria during the first wave (January to May) of the coronavirus disease 2019 (COVID-19) pandemic. We found that the implementation of individual travel history was essential for the estimation of SARS-CoV-2 movement via discrete phylogeography models. Without the additional information provided by the travel history metadata, the discrete trait analysis could not be fit to the data due to numerical instability. We also suggest that during the first wave of the COVID-19 pandemic in Australia, the primary driving force behind the spread of SARS-CoV-2 was viral importation from international locations. This case study demonstrates the necessity of robust genomic datasets supplemented with epidemiological metadata for generating accurate estimates from phylogeographical models in datasets that have significant sampling bias. For future work, we recommend the collection of metadata in conjunction with genomic data. Furthermore, we highlight the risk of applying phylogeographical models to biased datasets without incorporating appropriate metadata, especially when estimates influence public health policy decision making.

Phylogenomic investigation of an outbreak of fluoroquinolone-resistant Salmonella enterica subsp. enterica serovar Paratyphi A in Phnom Penh, Cambodia.

In early 2020, the Medical Biology Laboratory of the Pasteur Institute of Cambodia isolated an unusually high number of fluoroquinolone-resistant Salmonella enterica subspecies enterica serovar Paratyphi A strains during its routine bacteriological surveillance activities in Phnom Penh, Cambodia. A public-health investigation was supported by genome sequencing of these Paratyphi A strains to gain insights into the genetic diversity and population structure of a potential outbreak of fluoroquinolone-resistant paratyphoid fever. Comparative genomic and phylodynamic analyses revealed the 2020 strains were descended from a previously described 2013-2015 outbreak of Paratyphi A infections. Our analysis showed sub-lineage 2.3.1 had remained largely susceptible to fluoroquinolone drugs until 2015, but acquired chromosomal resistance to these drugs during six separate events between late 2012 and 2015. The emergence of fluoroquinolone resistance was rapidly followed by the replacement of the original susceptible Paratyphi A population, which led to a dramatic increase of fluoroquinolone-resistant blood-culture-confirmed cases in subsequent years (2016-2020). The rapid acquisition of resistance-conferring mutations in the Paratyphi A population over a 3 year period is suggestive of a strong selective pressure on that population, likely linked with fluoroquinolone use. In turn, emergence of fluoroquinolone resistance has led to increased use of extended-spectrum cephalosporins like ceftriaxone that are becoming the drug of choice for empirical treatment of paratyphoid fever in Cambodia.

An ISO-certified genomics workflow for identification and surveillance of antimicrobial resistance.

Realising the promise of genomics to revolutionise identification and surveillance of antimicrobial resistance (AMR) has been a long-standing challenge in clinical and public health microbiology. Here, we report the creation and validation of abritAMR, an ISO-certified bioinformatics platform for genomics-based bacterial AMR gene detection. The abritAMR platform utilises NCBI's AMRFinderPlus, as well as additional features that classify AMR determinants into antibiotic classes and provide customised reports. We validate abritAMR by comparing with PCR or reference genomes, representing 1500 different bacteria and 415 resistance alleles. In these analyses, abritAMR displays 99.9% accuracy, 97.9% sensitivity and 100% specificity. We also compared genomic predictions of phenotype for 864 Salmonella spp. against agar dilution results, showing 98.9% accuracy. The implementation of abritAMR in our institution has resulted in streamlined bioinformatics and reporting pathways, and has been readily updated and re-verified. The abritAMR tool and validation datasets are publicly available to assist laboratories everywhere harness the power of AMR genomics in professional practice.

Improved Genome Sequence of Australian Methicillin-Resistant Staphylococcus aureus Strain JKD6159.

Staphylococcus aureus strain JKD6159 represents a prominent community-acquired methicillin-resistant S. aureus (MRSA) clone in Australia. Here, we report an improved assembly of the original S. aureus JKD6159 genome sequence. By using deep sequencing with multiple technologies combined with carefully curated assembly and polishing, we believe the assembly to contain zero errors.

Genomic characterisation reveals a dominant lineage of SARS-CoV-2 in Papua New Guinea.

The coronavirus disease pandemic has highlighted the utility of pathogen genomics as a key part of comprehensive public health response to emerging infectious diseases threats, however, the ability to generate, analyse, and respond to pathogen genomic data varies around the world. Papua New Guinea (PNG), which has limited in-country capacity for genomics, has experienced significant outbreaks of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) with initial genomics data indicating a large proportion of cases were from lineages that are not well defined within the current nomenclature. Through a partnership between in-country public health agencies and academic organisations, industry, and a public health genomics reference laboratory in Australia a system for routine SARS-CoV-2 genomics from PNG was established. Here we aim to characterise and describe the genomics of PNG's second wave and examine the sudden expansion of a lineage that is not well defined but very prevalent in the Western Pacific region. We generated 1797 sequences from cases in PNG and performed phylogenetic and phylodynamic analyses to examine the outbreak and characterise the circulating lineages and clusters present. Our results reveal the rapid expansion of the B.1.466.2 and related lineages within PNG, from multiple introductions into the country. We also highlight the difficulties that unstable lineage assignment causes when using genomics to assist with rapid cluster definitions.

Niche-specific genome degradation and convergent evolution shaping Staphylococcus aureus adaptation during severe infections.

During severe infections, Staphylococcus aureus moves from its colonising sites to blood and tissues and is exposed to new selective pressures, thus, potentially driving adaptive evolution. Previous studies have shown the key role of the agr locus in S. aureus pathoadaptation; however, a more comprehensive characterisation of genetic signatures of bacterial adaptation may enable prediction of clinical outcomes and reveal new targets for treatment and prevention of these infections. Here, we measured adaptation using within-host evolution analysis of 2590 S. aureus genomes from 396 independent episodes of infection. By capturing a comprehensive repertoire of single nucleotide and structural genome variations, we found evidence of a distinctive evolutionary pattern within the infecting populations compared to colonising bacteria. These invasive strains had up to 20-fold enrichments for genome degradation signatures and displayed significantly convergent mutations in a distinctive set of genes, linked to antibiotic response and pathogenesis. In addition to agr-mediated adaptation, we identified non-canonical, genome-wide significant loci including sucA-sucB and stp1. The prevalence of adaptive changes increased with infection extent, emphasising the clinical significance of these signatures. These findings provide a high-resolution picture of the molecular changes when S. aureus transitions from colonisation to severe infection and may inform correlation of infection outcomes with adaptation signatures.

The bacterium Staphylococcus aureus lives harmlessly on our skin and noses. However, occasionally, it gets into our blood and internal organs, such as our bones and joints, where it causes severe, long-lasting infections that are difficult to treat. Over time, S. aureus acquire characteristics that help them to adapt to different locations, such as transitioning from the nose to the blood, and avoid being killed by antibiotics. Previous studies have identified changes, or 'mutations', in genes that are likely to play an important role in this evolutionary process. One of these genes, called accessory gene regulator (or agr for short), has been shown to control the mechanisms S. aureus use to infect cells and disseminate in the body. However, it is unclear if there are changes in other genes that also help S. aureus adapt to life inside the human body. To help resolve this mystery, Giulieri et al. collected 2,500 samples of S. aureus from almost 400 people. This included bacteria harmlessly living on the skin or in the nose, as well as strains that caused an infection. Gene sequencing revealed a small number of genes, referred to as 'adaptive genes', that often acquire mutations during infection. Of these, agr was the most commonly altered. However, mutations in less well-known genes were also identified: some of these genes are related to resistance to antibiotics, while others are involved in chemical processes that help the bacteria to process nutrients. Most mutations were caused by random errors being introduced in to the bacteria's genetic code which stopped genes from working. However, in some cases, genes were turned off by small fragments of DNA moving around and inserting themselves into different parts of the genome. This study highlights a group of genes that help S. aureus to thrive inside the body and cause severe and prolonged infections. If these results can be confirmed, it may help to guide which antibiotics are used to treat different infections. Furthermore, understanding which genes are important for infection could lead to new strategies for eliminating this dangerous bacterium.

State-wide genomic epidemiology investigations of COVID-19 in healthcare workers in 2020 Victoria, Australia: Qualitative thematic analysis to provide insights for future pandemic preparedness.

Background: COVID-19 has affected many healthcare workers (HCWs) globally. We performed state-wide SARS-CoV-2 genomic epidemiological investigations to identify HCW transmission dynamics and provide recommendations to optimise healthcare system preparedness for future outbreaks. Methods: Genome sequencing was attempted on all COVID-19 cases in Victoria, Australia. We combined genomic and epidemiologic data to investigate the source of HCW infections across multiple healthcare facilities (HCFs) in the state. Phylogenetic analysis and fine-scale hierarchical clustering were performed for the entire dataset including community and healthcare cases. Facilities provided standardised epidemiological data and putative transmission links. Findings: Between March-October 2020, approximately 1,240 HCW COVID-19 infection cases were identified; 765 are included here, requested for hospital investigations. Genomic sequencing was successful for 612 (80%) cases. Thirty-six investigations were undertaken across 12 HCFs. Genomic analysis revealed that multiple introductions of COVID-19 into facilities (31/36) were more common than single introductions (5/36). Major contributors to HCW acquisitions included mobility of staff and patients between wards and facilities, and characteristics and behaviours of patients that generated numerous secondary infections. Key limitations at the HCF level were identified. Interpretation: Genomic epidemiological analyses enhanced understanding of HCW infections, revealing unsuspected clusters and transmission networks. Combined analysis of all HCWs and patients in a HCF should be conducted, supported by high rates of sequencing coverage for all cases in the population. Established systems for integrated genomic epidemiological investigations in healthcare settings will improve HCW safety in future pandemics. Funding: The Victorian Government, the National Health and Medical Research Council Australia, and the Medical Research Future Fund.

SARS-CoV-2 infection results in immune responses in the respiratory tract and peripheral blood that suggest mechanisms of disease severity.

Respiratory tract infection with SARS-CoV-2 results in varying immunopathology underlying COVID-19. We examine cellular, humoral and cytokine responses covering 382 immune components in longitudinal blood and respiratory samples from hospitalized COVID-19 patients. SARS-CoV-2-specific IgM, IgG, IgA are detected in respiratory tract and blood, however, receptor-binding domain (RBD)-specific IgM and IgG seroconversion is enhanced in respiratory specimens. SARS-CoV-2 neutralization activity in respiratory samples correlates with RBD-specific IgM and IgG levels. Cytokines/chemokines vary between respiratory samples and plasma, indicating that inflammation should be assessed in respiratory specimens to understand immunopathology. IFN-α2 and IL-12p70 in endotracheal aspirate and neutralization in sputum negatively correlate with duration of hospital stay. Diverse immune subsets are detected in respiratory samples, dominated by neutrophils. Importantly, dexamethasone treatment does not affect humoral responses in blood of COVID-19 patients. Our study unveils differential immune responses between respiratory samples and blood, and shows how drug therapy affects immune responses during COVID-19.

Genomic Evidence of In-Flight SARS-CoV-2 Transmission, India to Australia, April 2021.

Epidemiologic and genomic investigation of SARS-CoV-2 infections associated with 2 repatriation flights from India to Australia in April 2021 indicated that 4 passengers transmitted SARS-CoV-2 to >11 other passengers. Results suggest transmission despite mandatory mask use and predeparture testing. For subsequent flights, predeparture quarantine and expanded predeparture testing were implemented.

Multi-site implementation of whole genome sequencing for hospital infection control: A prospective genomic epidemiological analysis.

Background: Current microbiological methods lack the resolution to accurately identify multidrug-resistant organism (MDRO) transmission, however, whole genome sequencing can identify highly-related patient isolates providing opportunities for precision infection control interventions. We investigated the feasibility and potential impact of a prospective multi-centre genomics workflow for hospital infection control. Methods: We conducted a prospective genomics implementation study across eight Australian hospitals over 15 months (2017,2018), collecting all clinical and screening isolates from inpatients with vanA VRE, MRSA, ESBL Escherichia coli (ESBL-Ec), or ESBL Klebsiella pneumoniae (ESBL-Kp). Genomic and epidemiologic data were integrated to assess MDRO transmission. Findings: In total, 2275 isolates were included from 1970 patients, predominantly ESBL-Ec (40·8%) followed by MRSA (35·6%), vanA VRE (15·2%), and ESBL-Kp (8·3%).Overall, hospital and genomic epidemiology showed 607 patients (30·8%) acquired their MDRO in hospital, including the majority of vanA VRE (266 patients, 86·4%), with lower proportions of ESBL-Ec (186 patients, 23·0%), ESBL-Kp (42 patients, 26·3%), and MRSA (113 patients, 16·3%). Complex patient movements meant the majority of MDRO transmissions would remain undetected without genomic data.The genomics implementation had major impacts, identifying unexpected MDRO transmissions prompting new infection control interventions, and contributing to vanA VRE becoming a notifiable condition. We identified barriers to implementation and recommend strategies for mitigation. Interpretation: Implementation of a multi-centre genomics-informed infection control workflow is feasible and identifies many unrecognised MDRO transmissions. This provides critical opportunities for interventions to improve patient safety in hospitals. Funding: Melbourne Genomics Health Alliance (supported by State Government of Victoria, Australia), and National Health and Medical Research Council (Australia).

Optimising genomic approaches for identifying vancomycin-resistant Enterococcus faecium transmission in healthcare settings.

Vancomycin-resistant Enterococcus faecium (VREfm) is a major nosocomial pathogen. Identifying VREfm transmission dynamics permits targeted interventions, and while genomics is increasingly being utilised, methods are not yet standardised or optimised for accuracy. We aimed to develop a standardized genomic method for identifying putative VREfm transmission links. Using comprehensive genomic and epidemiological data from a cohort of 308 VREfm infection or colonization cases, we compared multiple approaches for quantifying genetic relatedness. We showed that clustering by core genome multilocus sequence type (cgMLST) was more informative of population structure than traditional MLST. Pairwise genome comparisons using split k-mer analysis (SKA) provided the high-level resolution needed to infer patient-to-patient transmission. The more common mapping to a reference genome was not sufficiently discriminatory, defining more than three times more genomic transmission events than SKA (3729 compared to 1079 events). Here, we show a standardized genomic framework for inferring VREfm transmission that can be the basis for global deployment of VREfm genomics into routine outbreak detection and investigation.

Daptomycin Resistance Occurs Predominantly in vanA-Type Vancomycin-Resistant Enterococcus faecium in Australasia and Is Associated With Heterogeneous and Novel Mutations.

Healthcare associated infections caused by vancomycin-resistant Enterococcus faecium (VREfm) have a major impact on health outcomes. VREfm is difficult to treat because of intrinsic and acquired resistance to many clinically used antimicrobials, with daptomycin being one of the few last line therapeutic options for treating multidrug-resistant VREfm. The emergence of daptomycin-resistant VREfm is therefore of serious clinical concern. Despite this, the impact that daptomycin-resistant VREfm have on patient health outcomes is not clearly defined and knowledge on the mechanisms and genetic signatures linked with daptomycin resistance in VREfm remains incomplete. To address these knowledge gaps, phenotypic daptomycin susceptibility testing was undertaken on 324 E. faecium isolates from Australia and New Zealand. Approximately 15% of study isolates were phenotypically resistant to daptomycin. Whole genome sequencing revealed a strong association between vanA-VREfm and daptomycin resistance, with 95% of daptomycin-resistant study isolates harbouring vanA. Genomic analyses showed that daptomycin-resistant VREfm isolates were polyclonal and carried several previously characterised mutations in the liaR and liaS genes as well as several novel mutations within the rpoB, rpoC, and dltC genes. Overall, 70% of daptomycin-resistant study isolates were found to carry mutations within the liaR, rpoB, rpoC, or dltC genes. Finally, in a mouse model of VREfm bacteraemia, infection with the locally dominant daptomycin-resistant clone led to reduced daptomycin treatment efficacy in comparison to daptomycin-susceptible E. faecium. These findings have important implications for ongoing VREfm surveillance activities and the treatment of VREfm infections.